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Abstract 

We study general properties of attractors for tachyonic potential chameleon scalar-field model 
which possess cosmological scaling solutions. An analytic formulation is given to obtain fixed points 
with a discussion on their stability. The model predicts a dynamical equation of state parameter 
with phantom crossing behavior for an accelerating universe. We constrain the parameters of the 
model by best fitting with the recent data-sets from supernovae and simulated data points for 
redshift drift experiment generated by Monte Carlo simulations. 
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1. INTRODUCTION 



Recently, the observations of high redshift type la supernovae and the surveys of clusters 
of galaxies [l| reveal the universe accelerating expansion and that the density of matter 
is very much less than the critical density. Also the observations of Cosmic Microwave 
Background (CMB) anisotropies indicate that the universe is fiat and the total energy density 
is very close to the critical one j2|. While the above observational data properly complete 
each other, the dynamical dark energy proposal as an interesting possibility may arise to 
explain the observational constraints |3|. 

The two fine tuning problem and cosmic coincidence problem are the most serious issues 
with regards to the dark energy (DE) models and of the most frequently used approach to 



5 scenario by employing scalar field models 
J]. Tracker models are independent of initial 



moderate these problems is the tracker field D 
which exhibit scaling solutions, for example see 
conditions used for field evolution but do require the tuning of the slope of the scalar field 
potential. During the scaling regime, the scalar field energy density is of the same order of 
magnitude as the background energy density. The scaling solutions as dynamical attractors 
can considerably resolve the two above mentioned problems. Furthermore, by investigating 
the nature of scaling solutions, one can determine whether such behavior is stable or just a 
transient feature and explore the asymptotic behavior of the scalar field potential j^. 

In studying the canonical scalar field models with exponential potential there exist scaling 
attractor solutions ^|. In quintessence dark energy model, there are two scaling solutions. 
The first one is fluid-scalar field scaling solution, which remains subdominant for most of the 
cosmic evolution. It is necessary that the scalar field mimics the background energy density 
(radiation/ matter) in order to respect the nucleosynthesis constraint and can also solve the 
fine-tuning problem of initial conditions. The second one is scalar field dominated scaling 
solution, which is a late time attractor and gives rise to the accelerated expansion. Since 
the fluid-scalar field scaling solution is non-accelerating, we need an additional mechanism 
exit from the scaling regime so as to enter the scalar field dominated scaling solution at late 
times. For the discussion on^the exiting mechanism one can refer to [3]. 



In scalar-tensor theories 18 



17| , interaction of the scalar field with matter ( for example in 



chameleon cosmology) 18|- 20| and the presence of the tachyon potential in the formalism 



2l|- 22| separately are used to interpret the late time acceleration. From a chameleonic 
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point of view, since the field mimics tfie background radiation/matter field, subdominant 
for most of the evolution history except at late times when it becomes dominant, it may 
be regarded as a cosmological tracker field {23]. For tachyon dark energy J24, 27| as well 
as other scalar-tensor models such as chameleonic Brans-Dicke cosmology 28|, the scaling 
solutions have been investigated separately. 

Both Chameleon and tachyon fields have been motivated separately from modified gravity 



and string theory 



29|- 



30| and explain the quintessence cosmological scenario and the late 



time universe acceleration. However in order to interpret also the phantom behavior of dark 
energy, we integrate both models into one to investigate in addition to the scaling solutions 
and late time acceleration of the universe, the phantom crossing and attractor behavior of 
;he solutions which are also obtained in the context of non- minimally coupled scalar fields 



3l| and are consistent with the observations. 



The scalar field in our model plays the role of chameleon field coupled with the matter 
lagrangian and also by its presence in tachyonic potential can be regarded as tachyon field. 
The scalar field coupled with the matter lagrangian with a prefect fluid, pm = IPm, while 
the tachyonic scalar field candidates for dark energy. 

The well-known geometric variables, i.e. Hubble and deceleration parameters at the 
present time are used to explain the acceleration expansion of the universe. However, con- 
sidering the increased accuracy of the observational data during the last few years and 
generality of the DE models, new geometrical variables are introduced to differentiate these 
models and better fit the observational data. In this regard, the authors in 32[ proposed 
a cosmological diagnostic pair {s,r}, called statefinder, to differentiate the expansion dy- 
namics with higher derivatives of the scale factor and is a natural next step beyond the 
well known geometric variables. The statefinder pair has been used to explore a series of 
dark energy and cosmological models, including A cold dark matter (ACDM), quintessence, 
coupled quintessence, Chaplygin gas, holographic dark energy models, braneworld models, 



and so on 



33 



The Cosmological Redshift Drift (CRD) test which maps the expansion of the universe 
directly is also examined in here. We assume that the universe is homogeneous and isotropic 
at the cosmological scales 36| . The test is based on very simple and straightforward physics. 



Observationally, it is a very challenging task and requires technological breakthroughs, for 
more details see {s^I-ImI- we also verify our model by comparing the distance modulus mea- 
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surements versus redshift derived in the model with the data obtained from the observations 
of type la supernovae and utihzing the method. 

In the following, we study the detailed evolution of the model and the attractor property 
of its solution. In Sec. 2, we derive the field equations. In Sec. 3, we obtain the autonomous 
equations of the model and by using the phase plane analysis qualitatively investigate the 
dynamic of the system and the existence of a late time attractor solution. We also examine 
the behavior of the EoS parameter of the model and also perform a statefinder diagnostic for 
the model and analyze the evolving trajectories of the model in the statefinder parameter 
plane. Section 4 is devoted to perform cosmological tests. In sec. 5 we present summary 
and remarks. 



2. THE MODEL 



The chameleon gravity with a tachyonic potential is given by, 

S = J - V{<P)^1 - d,<Pd^cP + f{<P)Cm]V^dx\ (1) 

where R is Ricci scalar, G is the newtonian constant gravity, and the second term in the 
action is tachyon potential. Unlike the usual Einstein-Hilbert action in chameleon cosmology, 
the matter Lagrangian Cm is modified as f{(j))C„i , where /(</>) is an analytic function of the 
scalar field. The last term in Lagrangian brings about the nonminimal interaction between 
the matter and chameleon field. The variation of action ([T]) with respect to the metric tensor 
components in a spatially flat FRW cosmology yields the field equations: 

3//^ = P./ + ^i£L, (2) 



y 1- 

2H + ?,H^ = -^pmf + V{ct))^Jl - 02, (3) 

where we put SnG = c = h = 1 and H = ^ with a is the scale factor of the universe. We 
also assume a perfect fluid with p„i = jpm- Note that in here 7 is the EoS parameter for 
nonminimally coupled chameleon field with the matter field in the universe. Variation of 
the action ([1]) with respect to the scalar field (/) provides the wave equation for chameleon 
field as 

+ (1 - 02)(3//0 + ^) = - 02)1^^^ (4) 



where prime indicates differentiation with respect to (f) and e = From equations ([2])-([3]) 

one arrives at the conservation equation, 

(Pm/) + 'iH{l + 7)pm,/ = -epmf- (5) 

From equations (E]) and (E]), by defining the effective energy density and pressure, pe// 
and Peffi one can identify an effective EoS parameter as 



Peff Pmf + 



— Peff _ it^'"-' • \-r / y - -r 
^eff = - . . V(^) 



In the next section we study the stabihty of the system in the frame work of phase-space 
trajectory analysis. 

3. STABILITY ANALYSIS-PERTURBATION AND PHASE SPACE 

The structure of the dynamical system can be studied via phase plane analysis, by intro- 
ducing the following dimensionless variables, 

Pmf V ■ 1 

We consider that /(0) = /oexp(5i0) and V{(/)) = Voexp(520) where 6i and 62 are dimen- 
sionless constants characterizing the slope of potential V{(j)) and coupling field f{4>). The 
cosmological models with such exponential functions have been known lead to interesting 



physics in a variety of context, ranging from existence of accelerated expansions 42| to 



cosmological scaling solutions 43|-4^. In 



^articular, the exponential forms of /(^ and 



V{(p) are motivated by chameleon models 29|] and also from stability considerations j45j. In 



addition, attractor solutions with exponential functions may lead to cosmic acceleration for 
natural values of model parameters 461 ] . 



By using (l2])-([5]), the equations for the new dynamical variables are, 
dx 3yz^ 

— = y[3{l + 7)x + 3^== + — ], (9) 



„ „ .12. 

dN 





7T "T 
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(10) 



^ = ^[-(1 + j)x + - 1, (11) 



where = /n(a). By using the Fridmann constraint equation ([2]) in terms of the new 
dynamical variables: 



the equations (l8l)- f[TT|) reduce to, 

= x[3{l — x)z'^ — diewz], (13) 
^ = (^^-l)(3. + + — (14) 
1^ = + 7)1 + ^(1 (15) 

In stability formalism, by simultaneously solving = 0, = and to' = 0, where prime 
from now on means derivative with respect to /n(a), the fixed points (critical points) can be 
obtained. The critical points that may explicitly depend on the cosmological and stability 
parameters 7, 5i and 82 are illustrated in Table I. Substituting linear perturbations x' — )■ x' + 
z' — > z'+5z\ w' —J- w'+6w' about the critical points into the three independent equations 
(fT3|) - (fT5!) . to the first orders in the perturbations, gives us three eigenvalues Xi{i = 1,2,3) 
which has to be negative as a requirement by stability method. 

TABLE I: critical points 



points PI P2 


P3 


P4 


P5 


a; 


x{t) 


_i 

7 


_i 

7 


z 1 -1 





1 


-1 


w 





-12(1+7) 


12(1+7) 


<5i7(-l+37) 


<5i7(-l+37) 



In the following, the nature of the five critical points are given with the stability conditions 

as 
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PI : 


: Aipi 


= 6, A2P1 


= 3, A3P1 = 


P2 : 


: AiP2 


= 6, A2P2 


= 3, A3P2 = 


P3 : 


'■ AiP3 


= 0, A2P3 


= ~3, A3P3 



3 

-, unstable 
2' 

3 

-, unstable 
2' 

= ^x(t)(l + 



p.p., _ 3(1 + v/l + 2f + 27) _ -3(-l + + + 27) 

-r4 , i-'O : AiP4 5 — , A2P4,5 



A 



3P4,5 



27 ' ' 27 

6(2^17 + 3517^ - A62 - 4(527 - ^i) 
7(-l + 37)5i ' 



stable for(52 < 7(37 - -1< 7 < 0, ^ > 0, (16) 

4 70i(d7 — Ij 

or 52 > 7(37-1)^1, -1<7<0, ^ < 

4 7di(37-l) 

As can be seen, the critical points PI, P2 and P3 are unstable and the points P4 and 
P5 are stable for the given conditions on 61, 62 and 7. It is interesting to note that in order 
for the critical points p4 and P5 to be stable, in addition to the conditions on the stability 
parameters di and 62, the EoS parameter of the matter field has to be between zero and 
— 1. This means that the system is stable for the chameleon field interaction with dark 
energy. In the following we will numerically investigate the models with specific choice of 
EoS parameter of the matter in the universe. In Fig. 1, the attraction of the trajectories 
to the critical points, P4 and P5 in the 3-dim phase plane is shown for the given initial 
conditions. 
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Fig. 1: The attractor property of the dynamical system in the 3-dim phase plane 
for the case 7 = —0.6, 61 = —12 and 62 = —10 and different initial conditions 

Fig. 2 shows the projection of the 3-dim phase space of P4 and P5 onto the 2-dim phase 
plane: 
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stable • 




1 1.5 




stable • 





Fig. 2: The 2-dim phase plane corresponding to the critical points P4(top figures) and 
P5 (bottom figures). We have 7 = —0.6, 61 = —12 and 62 = —10 
I.C.s:top)left)blue:x(0) = 2,«;(0) = -4,red:x(0) = 1.5,u;(0) = 2,black:x(0) = 2,«;(0) = 0.9 
,green:x(0) = 2.8,w{0) = 0.7,cyan:x(0) = 2.5, w{0) = 0.06,top)right)blue: x(0) = 2, 
z{0) = 0.1,red:a;(0) = 1.5,^(0) = 0.2,black:x(0) = 2,^(0) = 0.6,green:x(0) = 2.8, 
^(0) = 0.6,cyan:x(0) = 2.5,^(0) = 0.9,bottom)left)blue:x(0) = 2,w{0) = -4,red:2;(0) = 2.5, 
w{0) = -3.5,black:2;(0) = 2,w(0) = -3,green:x(0) = 3.8, w{0) = -3.5,cyan:x(0) = 2, 
w{0) = -2,gray:x(0) = 2.5, w{0) = -3,bottom)right)blue: x(0) = 2, z{0) = 0.3, 
red:x(0) = 1.5, z{0) = 0.8,black:x(0) = 2,^(0) = -0.9,green: x(0) = 2.8, z{0) = 0.9, 
cyan: x(0) = 2.5, ^(0) = -l,gray: x(0) = 2.5, x(0) = -0.9, 

The particular choices of the model parameters 61, 62, 7 and initial conditions are just 
examples in expression (fT6|) within the range of stable critical points P4 and P5 (the at- 
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tractors). If we chose any other values for these parameters within the range, and initial 
conditions for the dynamical variables near the fixed points P4 and P5, we still get the same 
behavior for the trajectories. That is, with small perturbation, the trajectories approach 
the same stable critical points as shown in figures 1 and 2. 

Two of the cosmological parameters which relates the dynamics of the universe with the 
observational data are the EoS and deceleration parameters of the cosmological model. In 
terms of the new dynamical variables in our model they are given by. 



UJ. 



eff 



jx — yVl — z"^, (17) 



H _ 3 



q = -l- — = -l + -(7x - yVl -z^ + 1). (18) 

Moreover, the statefinder parameters {r, s}, the effective squared sound speed and the power 
p in the scale factor a a in terms of the dynamical variables are, 

r = ^z{2 - z')A + ^-^{1 - x)B + |x7(7 + 1) - ^(1 - ^)^" + 1> (19) 



lz{2 - z^)A + f^(l - x)i3 + |x7(7 + 1) - 1(1 - x)z^ 
1(72; - y\/l - z2) 



(20) 



9 '^'eff 'IX' + X' {1 — Z^) + 2(1 — x)zz' , X/. Ox 

3(a;e// + l) 3((1 - x)(l - z^) - 1 - 7a;) 

P = 1 , (22) 

3(7X — yv 1 — z"^ + 1) 

where A = 3(1 — x)z + S2Uj{1 — x) — Siexw and B = {1 + 'y)e6ixV^ — z"^ — 52z'^. 
For our model these parameters are presented in Table II, 



TABLE II: Properties of the critical points 



points 



^eff 



acceleration 



PI 

P2 

P3 x(7 + l) 

P4 -00 

P5 -00 



1/2 
1/2 

1 ix(7 + l) 
-1 
-1 






1 x(7 + l)-l 
-1 
-1 



1 
1 

1 + |x7(7 + 1) 



_Z _L 9^+7 
2 2V^P~ 

7 I 9/1+7 
2 2\~y^ 



7) 

7) - 36C 8C 



00 
00 

x(7+l)-l 
1+7 
7^ 



+ 7+1 



1+7 
7^ 



+ 7 + 1 



No 
No 

Conditional 
Yes 
Yes 
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where C is 



852(l+7)(l+i) 



(Si7(37-1) 

Note that with the disturbances propagate with the sound speed, the larger the sound 
speed, larger the range of influence. As long as the sound speed is finite, still the influence 
moves locally . For infinite sound speed, then disturbances can propagate instantaneously, 
lence the system becomes non-local and this is the case for incompressible fluids (5pe// = 0) 
47( 1 ■ As can be seen from Table II and Fig.3, the universe starts from an unstable state 
(critical point PI and P2) in the past when uo^ff = 0. It crosses the phantom divide line 
in near future at about 0.002 < z < 0.004 and eventually become tangent to the phantom 
divide line at far future when approaches the stable critical point P5. Note that for the 
stable critical points P4 and P5, the squared sound speed is minus infinity as expected from 
the behavior of Ueff (Fig. 3) approaching phantom divide line in far future. 

In our model, based on stability analysis, the critical points P4 and P5 are stable. That 
is, for a small perturbation of the parameters, the universe starts from the unstable state say 
PI and approaches the stable state P4 or P5. However, since the stable critical points are 
phantom state with Ueff = — 1, one would expect become minus infinity which is related 
to the nature of phantom energy. 

Also, from Fig. 3, the universe begins from an unstable state in the past with Ueff = 
and reaches a stable state in the future (phantom state) with Ueff = —1. Our result is in 
agreement with the work in 48| where from stability point of view some of the fixed points 
in the cosmological model are stable whereas they possess negative value for 49 1 



4r 




Ma) 



Fig. 3: The evolution of EoS parameter vs ln{a). ICs. x{0) = 1.6, z{0) = 1, w{0) = 2, 
e = 0.7 and 62 = -10, (black): 61 = -15, (red): di = -12, (blue): Si = -16. 



Fig. 4 shows the statefinder diagrams {s, q} and {r, q} evolutionary trajectories. From the 
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graph of the statefinder {r, q} all the trajectories for different values of stability parameter 
6i start from the standard cold dark energy (SCDM) in the past which is an unstable critical 
point and end their evolution to the critical point correspond to the state next to steady state 
(SS) in the future. The current values of the trajectories are shown and can be compared 
with the position of SCDM, LCDM and SS. From the graph of the statefinder {s, q} one 
observes that all the trajectories begin from the same unstable critical point and approach 
a state close to SS in the future. As can be seen, while the current values of the s and q are 
shown in the graph, all trajectories approach the stable critical point which is an attractor 
in the future. 

^1 • Today ~ SCDM 




-2-10 . . . . 

q q 



Fig. 4: Trajectories in the statefinder plane {r,q} and {s,q}. I.Cs. x{0) = 1.6, z{0) = 1, 
w{0) = 2, e = 0.7 and 62 = -10, (black): = -15, (red): 61 = -12, (blue): 61 = -16. 

In addition. Fig. 5 shows the statefinder diagram {r, s} evolutionary trajectory. From 
the graph we see that all the trajectories for different values of 61 commence evolving from 
different points in the past and continue towards the same stable critical point in the future 
near SS. 
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Fig. 5: Trajectories in the statefinder plane {r, s}, I.Cs. x(0) = 1.6, z(fi) = 1, w(0) = 2, 
e = 0.7 and §2 = -10, (black): 5i = -15, (red): Si = -12, (blue): 6i = -16. 

In Figs. 3,4 and 5, for effective EoS parameter and statefinders, we have selected stability 
parameters that are within the range for stable critical points P4 and P5. Specifically, we 
selected the same initial conditions and 6i for all of them and different 62, to show that 
as far as the chosen values are within the range, they all begin from the same points and 
approaches the attractors P4 or P5. If we would select other values for stability parameters 
within the range of stable critical point, still expecting to get the same behavior for the 
statefinder and effective EoS parameter. 

4. COSMOLOGICAL TEST 



We now examine our model with the observational data using the following cosmological 
tests |5oj-|53|. 



4.1. CPL, CRD and our model 



Following 



54l |. in Chevallier-Polarski-Linder (CPL) parametrization model one can use 



linearly approximated EoS parameter. 



Wcp/ - Wo 1 (« - 1) = Wo + Wi 

da 



1 + ^' 



(23) 
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where ooq is current value of the EoS and Ui = is its running factor. Using the above 

equation we can find the following equation for Hubble parameter, 



+ zf + (1 - + ^)3(i+-o+-i) X exp 



z 

-Swil-— 



(24) 



In CPL model the parametrization is fitted for different values of wq? and VLr 
On the other hand, the CRD can be extracted from 



z={l + z)H,-H{z), (25) 
which is known as Mc Vittie equation. This equation immediately leads to velocity drift 

„ = ci/„ - (26) 

1 + z 

By using equation ( 124|) in CPL model, the velocity drift with respect to the redshift can 
be obtained against observational data. In our model, from numerical computation one can 
obtain H{z) which can be used to evaluate i). To best fit the model for the parameters 
5i and 82 and the initial conditions x(0), -2(0), w(0), we use the simulated data points for 



redshift drift experiment generated by Monte Carlo simulations 53|] by employing the 
method. We constrain the parameters including the initial conditions by minimizing the 
function given as 



where the sum is over the three sets of data (8 points) for redshift drift experiments 53 |. 
In relation ( !27j) . vf^^ and vf'^ are the velocity drift parameters obtained from our model 
and from observation, respectively, and a is the estimated error of the vf"^. In our model 
the best fit values occur at 5i = -10.5, 62 = 3.5, x(0) = -0.6, z{0) = -0.2, w{0) = 1.09 
with Xmin — 1-27718546. In Fig. 6, the velocity drift, v, in our model is compared with the 
observational data for the obtained parameters and initial conditions using method in 
comparison with the CPL model and also ACDM model. As can be seen in our model for 
the best fitted stability parameters the model is in better agreement with the observational 
data in comparison with the CPL and ACDM models. 
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C0q= -2, cOj = 3.5 
C0p=-2.1,a)| =3.2 
ACDM model 



0.5 - _ . _ . . best — fitted model 




T ' 1 ' 1 ' 1 ' 1 ' T 

1 2 3 4 5 

z 

Fig. 6: The graph of velocity drift plotted as function of redshift for CPL model, 
our best-fitted model and ACDM . The stability parameter in our model: b\ = —10.5, 
62 = 3.5. ICs. x(0) = -0.6, ;^(0) = -0.2, w{0) = 1.09 and e = 0.7. 



4.2. The difference in the distance modulus, fi{z) 

The difference between the absolute and apparent luminosity of a distance object is given 
by, fi{z) = 25 + dlogiQdiiz) where the Luminosity distance quantity, dL^z) is given by 

In our model, from numerical computation one can obtain H{z) which can be used to 
evaluate yu(z). To best fit the model for the parameters 61 and 62 and the initial conditions 
x(0), z{0), w{0) with the most recent observational data, the Type la supernovea (SNe 
la), we employe the method. We constrain the model parameters including the initial 
conditions by minimizing the function given as 

XsNe{^U^2,x{0},Z[0),w{0)) = > 2 , (29) 

1=1 ' 

where the sum is over the SNe la sample. In relation ( 129|) . /i*'^^ and ^1°^"^ are the distance 
modulus parameters obtained from our model and from observation, respectively, and a is 
the estimated error of the /j,"^^. 
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For Gaussian distributed measurements, the likelihood function is proportional to 1"^. 
In Fig. 7 we show the 2-dim likelihood distributions for the parameters and ^2 as obtained 
from SNe la data. The parameters varies as 1.1 < b\ < 3.1 and —4 < ^2 < —5.6 with the 
maximum likelihood at b\ = 2.1, 62 = —4.8. 




Fig. 7: The graph of the 2-dim likelihood distribution for parameters 5i and 62- 
ICs. x(0) = 0.6,2(0) = 0.8,w{0) = l.l,dL{0) = 0. 

In Fig. 8, the distance modulus, fi{z), in our model is compared with the observational 
data for the obtained parameters and initial conditions using method. 




0.2 0.4 0.6 0.8 



Fig. 8: The graph of distance modulus fi{z) plotted as function of redshift, 

for 5i = 2.1, 52 = -4.8, e = 0.7. ICs. x(0) = 0.6,^(0) = 0.8,«;(0) = 1.1,^^(0) = 0. 
In our model the best-fitted model parameters occur at di = 2.1, 62 = —4.8, x(0) 



0.6, 
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^(0) = 0.8, w{0) = 1.1 with xLn = 541.6569729. 

5. SUMMARY AND CONCLUSION 

This paper is designed to study the attractor solutions of tachyonic potential chameleon 
cosmology by stability analysis and making use of the 3-dim phase space of the theory. The 
model characterized by the scalar field (j) , the scalar potential V{(f)), and the scalar function 
f{(p) nonminimally coupled to the matter lagrangian in the model. The phase space of 
the model is investigated using the scalar function f{(p) and potential V{(j)) in exponential 
form. The stability analysis gives the corresponding conditions for tracking attractor and 
determines the type of the universe behavior in the past and future. It has been shown that 
there are two stable critical points in the 3-dim phase space for the stability parameter 6i 
and 62- Together with the projection of the 3-dim phase space into 2-dim plane spaces, the 
trajectories for different stability parameters and I.Cs. are shown in Figs. 1 and 2. 

We then study the cosmological parameters such as effective EoS parameter, u^ff and 
statefinder parameters for the model in terms of the dynamical variables introduced in the 
stability section. It shows that the EoS parameter, We//, crosses the cosmological divide line 
in near future and become tangent to it in a stable state in far future. From Fig. 3, the 
universe starts from unstable state in the past with Ueff = and finally tends to the stable 
state in far future with Ueff = — 1. It is notable that for the stable critical points P4 and 
P5, the squared sound speed, cj. becomes minus infinity as Ueff tends to -1 in these stable 
states. The violation of the sound speed requirement to be positive is due to the property of 
the phantom states. From the statefinder graphs we see that in all scenarios with different 
stability parameters and initial conditions the universe begins from the standard cold dark 
energy (SCDM) unstable state and eventually tends to the stable state near SS state. It also 
shows that the universe with negative deceleration parameter is currently accelerating. To 
test our models, two cosmological tests, CRD, and the distance modulus are performed 
to compare our model with the CPL parametrization model, standard ACDM model and 
also with the observational data. In both tests, by employing the method we obtain a 
satisfactory fit to the observational data. Comparing the best-fitted velocity drift in our 
model with CPL and ACDM models shows that our model for the redshift around 2 and 
3.5 better matches the observational data. The best-fitted distance modulus ij,{z) with the 
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observational data, is also shown in Fig.8. 
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